Controllability of the discrete-spectrum Schrodinger 
equation driven by an external field 



00 

o 

O 

(N 



July 9, 2008 



Thomas Chambrion 

Institut Elic Cartan, UMR 7502 Nancy-Univcrsitc/CNRS/INRIA, BP 239, 54506 Vandccuvre-lcs-Nancy, France 
^\ ' Thomas.Chambrion@iecn.u-nancy.fr 

Paolo Mason 

lAC, CNR, Vialc Del Policlinico, 137 00161 Rome, Italy and 

Institut Elie Cartan, UMR 7502 INRIA/Nancy-Universitc/CNRS, BP 239, 54506 Vandceuvre-les-Nancy, France 

. p.mason@iac.cnr.it 
> . 

S2 ' Mario Sigalotti 



u 
o 



Institut Elie Cartan, UMR 7502 INRIA/Nancy-Universite/CNRS, BP 239, 54506 Vandceuvre-les-Nancy, France 

mario . sigalottiSinria. f r 

m ■ 

> ■ Ugo Boscain^ 

m ■ 

; Le2i, CNRS, Universite de Bourgogne, BP 47870, 21078 Dijon Cedex, France 
00 ' ugo.boscain@u-bourgogne.fr 

Abstract We prove approximate controllability of the bilinear Schrodinger 

O ' equation in the case in which the uncontrolled Hamiltonian has discrete non- 

00 

^ I resonant spectrum. The results that are obtained apply both to bounded or 

unbounded domains and to the case in which the control potential is bounded 
or unbounded. The method relies on finite-dimensional techniques applied to the 
Galerkin approximations and permits, in addition, to get some controllability prop- 
el I erties for the density matrix. Two examples are presented: the harmonic oscillator 
and the 3D well of potential, both controlled by suitable potentials. 
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Resume Nous montrons la controlabilite approchee de I'equation de Schrodinger 
bilineaire dans le cas oii I'hamiltonien non controle a un spectre discret et non- 
resonnant. Les resultats obtenus sont valables que le domaine soit borne ou non, et 
que le potentiel de controle soit borne ou non. La preuve repose sur des methodes 
de dimension finie appliquees aux approximations de Galerkyn du systeme. Ces 
methodes permettent en plus d'obtenir des resultats de controlabilite des matrices 
de densite. Deux exemples sont presentes, I'oscillateur harmonique et le puit de 
potentiel en dimension trois, munis de potentiels de controle adequats. 

The last author was partially supported by a FABER grant by C'onseil regional de Bourgogne 
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1 Introduction 



In this paper we study the controUabihty of the bihnear Schrodinger equation. Its importance 
is due to apphcations to modern technologies such as Nuclear Magnetic Resonance, laser 
spectroscopy, and quantum information science (see for instance [22, 29, 31, 38]). 

Many controllability results are available when the state space is finite dimensional, e.g., 
for spin systems or for molecular dynamics when one neglects interactions with highly excited 
levels (see for instance [8, 19]). When the state space is infinite-dimensional the controllability 
problem appears to be much more intricate. Some results are available when the control is 
the value of the wave function on some portion of the boundary or in some internal region of 
the domain (see [41] and references therein and the recent paper [39]). 

However, from the point of view of applications the case in which the control appears 
in the Hamiltonian as an external field is much more interesting, since the wave function is 
not directly accessible in experiments and because of the postulate of collapse of the wave 
function. For instance, in nuclear magnetic resonance the control is a magnetic field, in laser 
spectroscopy and in many applications of photochemistry the control is a laser or a source of 
hght. 

In this paper we consider the controllability problem for the following bilinear system 
representing the Schrodinger equation driven by one external field 



Here the wave function t/j evolves in an infinite-dimensional Hilbert space, Hq is a self-adjoint 
operator called drift Hamiltonian (i.e. the Hamiltonian responsible for the evolution when the 
external field is not active), u{t) is a scalar control function, and Hi is a self-adjoint operator 
describing the interrelation between the system and the external field. 

The reference case is the one in which the Hilbert space is L^{Q) where Q is either R'^ or 
a bounded domain of R"^, and equation (1.1) reads 



where A is the Laplacian (with Dirichlet boundary condition in the case in which Q is bounded) 
and V and W are suitably regular functions defined on Q. However the setting of the paper 
covers more general cases (for instance Q can be a Riemannian manifold and A the corre- 
sponding Laplace-Beltrami operator). Let us stress that the proposed approach allows to 
handle both cases where the control potential (i.e. Hi in (1.1) or W in (1.2)) is bounded or 
unbounded. Notice that in many situations the control potential happens to be unbounded. 
For instance if f2 = R*^ and the controlled external force depends on time, but is constant in 
space, then W is linear and hence unbounded. 

Besides the fact that one cannot expect exact controllability on the whole Hilbert sphere 
(see [10, 40]) and some negative result (in particular [28, 36]) only few approximate control- 
lability results are available and concern mainly special situations. It should be mentioned, 
however, that several results on efficient steering of the Schrodinger equation without any 
controllability assumptions are available, e.g. [11, 14, 23]. (For optimal control results for 
finite dimensional quantum systems see, for instance, [15, 16, 17, 26].) 



di/j 
I — 
dt 



it) = {Ho + u{t)Hi)m- 






2 



In [12, 13] Beauchard and Coron study the controllability of a quantum particle in a 
ID potential well with W{x) = x. Their results are highly nontrivial and are based on 
Coron's return method (see [18]) and Nash-Moser's theorem. In particular, they prove that 
the system is exactly controllable in the unit sphere of the Sobolev space if^ (implying in 
particular approximate controllability in L^). One of the most interesting corollaries of this 
result is exact controllability between eigenstates. 

A different result is given in [1], where adiabatic methods are used to prove approximate 
controllability for systems having conical eigenvalue crossings in the space of controls. 

Another controllability result has been proved by Mirrahimi in [27] using Strichartz esti- 
mates and concerns approximate controllability for a certain class of systems such that Q, = R'' 
and whose drift Hamiltonian has mixed spectrum (discrete and continuous). 

The aim of the present paper is to prove a general approximate controllability result for 
a large class of systems for which the drift Hamiltonian Hq has discrete spectrum. Our main 
assumptions are that the spectrum of Hq satisfies a non-resonance condition and that Hi 
couples each pair of distinct eigenstates of Hq. Such assumptions happen to be generic in a 
suitable sense, as it will be discussed in a forthcoming paper. 

We then apply the approximate controllability result to two classical examples, namely 
the harmonic oscillator and the 3D potential well, for suitable controlled potentials. 

Our method is new in the framework of quantum control and relies on finite-dimensional 
techniques applied to the Galerkin approximations. A difficult point is to deduce properties 
of the original infinite-dimensional system from its finite-dimensional approximations. For the 
Navier-Stokes equations this program was successfully conducted by Agrachev and Sarychev 
in the seminal paper [5] (see also [3, 35]). 

A key ingredient of the proof is a time reparametrization that inverts the roles of Hq 
and Hi as drift and control operator. This operation is crucial since it permits to exploit 
for the Galerkin approximation the techniques developed in [2] for finite-dimensional systems 
on compact semisimple Lie groups. The passage from the controllability properties of the 
Galerkin approximations to those of the infinite-dimensional system heavily relies on the fact 
that the dynamics preserve the Hilbert sphere. 

A feature of our method is that the infinite-dimensional system inherits, in a suitable sense, 
controllability results for the group of unitary transformations from those of the Galerkin 
approximations. This permits to extract controllability properties for the density matrix. Let 
us stress that, as it happens in finite dimension, controllability properties for the density 
matrix cannot in general be deduced from those of the wave function (see for instance [7]). 

The paper is organized as follows. In Section 2 we present the general functional analysis 
setting and we state our main result (Theorem 2.4) for the control system (1.1). In Section 
3 we show how this result applies to the Schrodinger equation (1.2) when Q is both bounded 
or unbounded. Section 4 contains the proof of Theorem 2.4 and an estimate of the minimum 
time for approximately steering the system between two given states, that holds even if the 
system itself is not approximately controllable. In Section 5 we extend Theorem 2.4 to the 
controlled evolution of the density matrix (Theorem 5.2). Finally in Section 6 we show how 
Theorem 2.4 and Theorem 5.2 can be applied to specific cases. In particular, we show how to 
get controllability results even in cases in which V does not satisfy the required non-resonance 
hypothesis, using perturbation arguments. 
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2 Mathematical framework and statement of the main 



result 



Hereafter N denotes the set of strictly positive integers. Definition 2.1 below provides the 
abstract mathematical framework that will be used to formulate and prove the controllability 
results later applied to the Schrodinger equation (1.2). The hypotheses under which (1.2) fits 
the abstract framework are discussed in Section 3. 

Definition 2.1. Let be a complex Hilbert space and f/ be a subset of R. Let A^B be 
two, possibly unbounded, operators on ^ with values in and denote by Di^A) and D{B) 
their domains. The control system {A, B, U) is the formal controlled equation 



We say that (A, 5, U) is a skew-adjoint discrete-spectrum control system if the following 
conditions are satisfied: (HI) A and B are skew-adjoint, (H2) there exists an orthonormal 
basis (0n)neN of =^ made of eigenvectors of A, (H3) 0„ G D{B) for every n G N. 

In order to give a meaning to the evolution equation (2.1), at least when u is constant, we 
should ensure that the sum A + uB is well defined. The standard notion of sum of operators 
seen as quadratic forms (see [20]) is not always applicable under the sole hypotheses (HI), 
(H2), (H3). An adapted definition of A + uB can nevertheless be given as follows: hypothesis 
(H3) guarantees that the sum A + uB is well defined on V = span{0„ | n G N}. Any skew- 
Hermitian operator C : V ^ M' admits a unique skew-adjoint extension ^(C). We identify 
A^uB with E{A\y + uB\y). 

Let us notice that when A + uB is well defined as sum of quadratic forms and is skew- 
adjoint then the two definitions of sum coincide. This happens in particular for the Schrodinger 
equation (1.2) in most physically significant situations (see Section 3). 

A crucial consequence of what precedes is that for every m G f/ the skew-adjoint operator 
A + uB generates a group of unitary transformations ^'^^^'^^^ jif j^f. In particular, the 
unit sphere S of satisfies e**-^"'"*'^^ (5) = S for every u E U and every t >0. 

Due to the dependence of the domain D{A + uB) on u, the solutions of (2.1) cannot in 
general be defined in classical (strong, mild or weak) sense. Let us mention that, in some 
relevant cases in which the spectrum of A has a nontrivial continuous component the solution 
can be defined as in [30, 34] by means of Strichartz estimates. 

We will say that the solution of (2.1) with initial condition ipo G and corresponding to 
the piecewise constant control u : [0,T] ^ U is the curve t h-> tpit) defined by 



(t) = Aij{t) +u{t)B^{t), 



u{t) G U. 



(2.1) 



o 



o ■ • • o e 



(2.2) 



where ^f^^ t« < t < J2i=i'^i '^i'^) — J2i=i'ti ^ < Si=i^«- Notice that such a ip^-) 
satisfies, for every n G N and almost every t G [0,T], the differential equation 



d 
dt 



(^(t),0n) 



(7A(t),(A + M(t)5)0„). 



(2.3) 
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Remark 2.2. The notion of solution introduced above makes sense in very degenerate situations 
and can be enhanced when B is bounded. Indeed, well-known results assert that in this case 
if M G L^{[0,T],U) then there exists a unique weak (and mild) solution tp G C{[0,T], J^) 
which coincides with the curve (2.2) when u is piecewise constant. Moreover, if ipQ G D{A) 
and u G C^([0,T], U) then ip is differentiable and it is a strong solution of (2.1). (See [10] and 
references therein.) 

Definition 2.3. Let {A, B, U) be a skew-adjoint discrete-spectrum control system. We say 
that [A, B, U) is approximately controllable if for every ipQ, ipi E S and every e > there exist 
A; G N, ti, . . . , tfc > and Ui, . . . , G U such that 

llVi - e**(^+"^^) o . . . o e*i(^+"^^)(^Ao)|| < e. 

Let, for every r?, G N, zA„ denote the eigenvalue of A corresponding to (A„ G R). The 
main result of the paper is the following. 

Theorem 2.4. Let 5 > and (A, i?, (0, 5)) he a skew-adjoint discrete- spectrum control system. 
If the elements of the sequence (A„+i — A„)„gN o^re Q-linearly independent and if {Bcpn, (pn+i) 7^ 
for every n G N, then [A, B, (0, 6)) is approximately controllable. 

Recall that the elements of the sequence (A„+i — A„)„gN are said to be Q-linearly inde- 
pendent if for every N eN and (gi, . . . , qjy) G \ {0} one has X^^^Li Qni^n+i — A„) 7^ 0. 

The condition {Bcpn, 4>n+i) 7^ 0, preferred here for the easiness of its expression, can be 
replaced by a weaker one (namely, (4.4)), as detailed in Remark 4.2. 

3 Discrete-spectrum Schrodinger operators 

The aim of this section is to recall some classical results on Schrodinger operators. In partic- 
ular we list here, among the numerous situations studied in the literature, some well-known 
sufficient conditions guaranteeing that the controlled Schrodinger equation (1.2) satisfies the 
assumptions of Definition 2.1. 

Theorem 3.1 ([21, Theorem 1.2.2]). Let Q be an open and bounded subset of Hf^ and V G 
L°°(r2,R). Then —A -|- V, with Dirichlet boundary conditions, is a self-adjoint operator with 
compact resolvent. In particular —A -\- V has discrete spectrum and admits a family of eigen- 
functions in H^{Q, R) fl IIq{Q, R) which forms an orthonormal basis of L'^{Q, C). 

Theorem 3.2 ([32, Theorems XIII.69 and XIII.70]). Let Q = R'^ and V e Ll^(R'^,R) be 
bounded from below and such that 

lim V{x) = -\-oo. 

|x|^oo 

Then —A -|- V , defined as a sum of quadratic forms, is a self-adjoint operator with compact 
resolvent. In particular —A-^V has discrete spectrum and admits a family of eigenf unctions in 
if^(R'^, R) which forms an orthonormal basis o/L^(R°', C). Moreover, for every eigenfunction 
(j) of —A -|- V and for every a > 0, x 1— e'*"^'"0(x) belongs to /.^(R'^, C). 
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In the following, we call controlled Schrodinger equation the partial differential equation 

i^{t, x) = (-A + V + uW)ip{t, x) 



where ip : I x Q ^ C, Q is an open subset of R'^, / is a subinterval of R and, in the case in 
which Q is bounded, ip\ixdn = 0. The correct functional analysis framework for this equation 
is specified below. 

The following corollary, which is a straightforward consequence of the results recalled 
above, states that the assumptions of Definition 2.1 are fulfilled by the operators appearing 
in the controlled Schrodinger equation under natural hypotheses. 

Corollary 3.3. Let Q be an open subset o/R*^, V, W be two real-valued functions defined on Q, 
and U be a subset o/R. Assume either that (i) Q is bounded, V, W belong to L°°{fl, R) or that 
(a) Q = IV^, V,W belong to //[^^(R'^, R), the growth of W at infinity is at most exponential 
and, for every u G U, lim||2.|j^+oo(^(a^) +uW{x)) = +oo and mi^^iid(y{x) +uW{x)) > — oo. 
Let be equal to L'^{Q,C) and D{A) be equal to H'^{fl,C) fl Hq{Q,C) in case (i) and to 
H'^{Q,C) in case (ii). Let, moreover, A be the differential operator — «(— A + V) and B be 
the multiplication operator —iW. Then {A,B,U) is a skew-adjoint discrete- spectrum control 
system, called the controlled Schrodinger equation associated with Q, V, W and U. 

Since the controlled Schrodinger equation is a skew-adjoint discrete-spectrum control sys- 
tem, it makes sense to apply Theorem 2.4 to it. The result is the following theorem. 

Theorem 3.4. Let Q, V, W and U satisfy one of the hypotheses (i) or (ii) of Corollary 3.3. 
Denote by {\k)ke'N the sequence of eigenvalues of —A -|- V and by {4>k)k£'N an orthonormal 
basis o/L^(f2, C) of corresponding real-valued eigenfunctions. Assume, in addition to (i) or 
(ii), that U contains the interval (0, 5) for some 5 > 0, that the elements of {\k+i — \k)k£'N are 
Q-linearly independent, and that W{x)(f)k4'k+i dx for every A; G N. Then the controlled 
Schrodinger equation associated with Q, V, W and U is approximately controllable. 

As stressed just after the statement of Theorem 2.4, the condition W{x)(l)k4>k+i dx 
could be replaced by a weaker one (see Remark 4.2). 

4 Proof of Theorem 2.4 

The proof of Theorem 2.4 is split in several steps. First, in Section 4.1 the controllability 
problem is transformed, thanks to a time-reparameterization, into an equivalent one where A 
and B play the role of controlled dynamics and drift, respectively. Then, in Section 4.2, we 
prove a controllability result for the Galerkin approximations of this equivalent system. In 
Section 4.3 we show how to lift the controllability properties from a Galerkin approximation 
to an higher-dimensional one. Section 4.4 makes the link between finite-dimensional and 
infinite-dimensional controllability properties and completes the proof. 

Finally, in Section 4.5, as a byproduct of the arguments of the proof, we get a lower bound 
on the minimum steering time. 
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4.1 Time-reparameterization 

First remark that, if u 7^ 0, = gt«((i/")A+s)_ theorem 2.4 is therefore equivalent to the 

following property: if the elements of the sequence (A„+i — A„)„gN are Q-linearly independent 
and if (-B0„,0„,+i) 7^ for every n G N, then for every 5, e > and every ipo,ipi G S there 
exist /c G N, ti, . . . , tfe > and Ui, . . . ,Uk > 6 such that 

In other words, the system for which the roles of A and B as drift and controlled field are 
inverted, namely, 

^{t) = + Bm, e u, (4.2) 

is approximately controllable provided that the control set U contains a half-line. The notion 
of solution of (4.2) corresponding to a piecewise constant control function is defined as in 
(2.2). 



4.2 Controllability of the Galerkin approximations 

Let, for every j,k G N, bjk = {B(f)j,(j)k) and ajk = {Acpj.cj)^) = iXj6jk- (Recall that \j G R 
and \ j G N} is the spectrum of A.) Define, for every n G N, the two complex- valued 
n X n matrices 

io,jk)i<j,k<n and 5'"' — {pjk)i<j,k<n- The Galerkyn approximation of 
(4.2) at order n (with respect to the basis (</'fc)fceN) is the finite-dimensional control system 

^ = uA^^'h + B^^'h, xeSn, u> 5, (S„) 
at 

where iS„ denotes the unit sphere of C". Notice that the system is well defined since, by 
construction, v4(") and 5^ are skew-Hermitian matrices. 

We say the (S„) is controllable if for every xo,xi G Sn there exist G N, ti, . . . , > 
and Ui, . . . ,Uk > 6 such that 

We recall that a n x n matrix C = {cjk)i<j^k<n is said to be connected if for every pair 
of indices j,k G {!,..., n} there exists a finite sequence ri, . . . , r; G {1, . . . , n} such that 
CjriCrir2 ' ' ' Cr, _ i r-j Cr, fc 7^ 0. (In the literature connected matrices are sometimes called einfach, 
or irreducible, or inseparable.) The following proposition is in the spirit of the controllability 
results obtained in [2] and [40]. 

Proposition 4.1. Let 21 = {<yjk)Jk=i! ^ ~ {f3jk)]k=i ^'"'^ skew- symmetric n x n matrices 
and assume that 21 is diagonal and *B is connected. Assume moreover that \ajj — akk\ 7^ 
|««/ — «mm| if {j, k} 7^ {/, m}. Then the control system (S) : x = u'Qlx + *Bx is controllable in 
Sn with piecewise constant controls u : H. ^ U, provided that U contains at least two points. 

Proof. For every 1 < j, k < n let cjk be the n x n matrix whose entries are all equal to zero 
except the one at line j and column k which is equal to 1. 
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Define for every p G N tfie iterated matrices commutator Mp = ad^(Q5). (Recall the usual 
notation adx(^) = [X, Y] = XY — YX for the adjoint operator associated with X.) A simple 
induction on p shows that the matrix Mp has the expression 

n 

l,m=l 

Fix two indices j ^ k such that I < j,k < n and Pjk 7^ 0. Since, by hypothesis, 
(-^kk Y 7^ — ftmrn)^ as soon as {j,k} ^ {/,?Ti}, there exists some polynomial P^^ 
with real coefficients such that Pjk[{c(jj — akkY) = 1 and Pjk[{aii — ammY) = for all 
{/, m} 7^ {j, k}, 1 < l,m < n. 

Let us define {ch)h as the coefficients of Pjk, i.e., Pjk{X) = J2h=o^hX^. Define moreover 
the matrix Njk = Ylh=o'^hM2h- By construction Njk = Y.lm.=i f^irnehnPjk{iau - a,nm)^) = 
PjkGjk — Pjk^kj- Therefore, the commutator [21, A^^-fc] is equal to {ajj — akk){Pjk€'jk + Pjk^kj) 
and so the Lie algebra generated by 21 and *B contains the two elementary anti-Hermitian 
matrices Ej^ = Cjk — e^j and Fj^ = i{ejk + e^j). 

Notice now that, for every 1 < j, k,h,m < n, ejkehm = ^kh^jm and therefore 

[Ejk.Ejk] = ^ii^Cjj — Ckk)- 

It follows from the definition of connected matrix and the relation [21, Ej^] = i{akk — (^jj)Pjk 
that the Lie algebra generated by 21 and 03 contains the matrices Ej^, Fj^ and i{ejj — e^k) for 
every j ^ k. Therefore 

su(n) C Lie(2t,?B). (4.3) 
Fix X G iS„ and consider the submersion V : SU{n) iS„ defined by V{g) = gx. Since 

Tv{g){Sn) = V,{TgSU{n)) = V,{m{n)g) = su{n)gx = su{n)V{g), 

then the evaluation at x = V{g) of the Lie algebra generated by 21 and ^ contains the whole 
space TxSn- Since for any u E U and any t G R the flow e*("2i+<8) ■ _^ jg volume- 
preserving then (E) is controllable (see [4, Cor. 8.6, Prop. 8.14, Th. 8.15]). □ 

The condition 7^ for every j G N appearing in Theorem 2.4 clearly ensures that 

every matrix 5^"^ is connected. Proposition 4.1, applied to A^"^ = 21 and -B'-"^ = 03, implies 
therefore that (S„) is controllable. 

Remark 4.2. In the following we can replace the assumption that fcj.j+i 7^ for every j G N 
with the weaker one that i?*-"-' is frequently connected, that is, 

Vj e'N,3k>j \ B^''^ is connected. (4.4) 

Notice, as a partial counterpart, that if there exists a nonempty and proper subset S of 
N such that for every j G S and G N \ S the coefficient bjj^ is equal to zero (i.e., the 
infinite-dimensional matrix (6/m)i,meN is non-connected) then the control system (2.1) is not 
approximately controllable. Indeed, the subspace span{0fc | G 2} is invariant for the 
dynamics of A + uB for every u E U and has nontrivial (invariant) orthogonal. 
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4.3 Approximate controllability in higher-dimensional projections 

Fix 6,e > and ipo, V^i G S. For every n G N, let n„ : be the orthogonal projection 

on the space span(0i, . . . , 0„) and n„ : C" be the map that associates to an element 

of the vector of its first n coordinates with respect to the basis (0m)meN- Choose n such 
that 

\\ilj,-Ii^{^ljj)\\<e for J = 0,1. (4.5) 

Thanks to (4.4) we can assume, without loss of generality, that (S„) is controllable. Let 
u : [0,T] — > (5, oo) be the piecewise constant control driving ^o/||^o|| to ^i/||^i|| where = 
Tiniil^j) for J = 0, 1. 

Let /i > be a constant which will be chosen later small enough, depending on T, n, 
and e. Notice that for every j G N the hypothesis that (pj belongs to D{B) implies that the 
sequence {bjk)k£'N is in l"^. It is therefore possible to choose N > n such that 

^ < yu, for every j = 1, . . . , n. (4.6) 

k>N 

If t I— > X(t) is a solution of (Sa?) corresponding to a control function U{-), then t i— > 
e-vm(''^X{t) = Y{t), where V{t) = fl^U{T)dT, is a solution of 



y(t) = e-^W^'"'i?(^)e^W^'"V(t). (6 



Let us represent the matrix e ^'(*)-4*^'^{A^)gi'(t)^*^' ^ where f(t) = j^u{r)dT^ in block form 
as follows 

^-.(t)A(-)5W^.(t)A(-) _ ^ M(".")(t) M("'^-")(t) \ 

e ^ ^ ~ \ M(^-"'")(t) M(^-"'^-")(t) )' ^ ' 

where the superscripts indicate the dimensions of each block. 

Claim 4.3. There exists a sequence of piecewise constant control functions Uk : [0, T] — ^ (5, oo) 
such that the sequence of matrix-valued curves 

where Vk{t) = j^Uk{T)dT, converges to 

t^M(t)-( ^^^'^^^'^ 



in the following integral sense 



t pt 
Mk{T)dT^ I M{T)dT {Al 



'0 ^0 

as k ^ oo uniformly with respect to t & [0,7"]. 

Proof. We will prove the claim taking v{-) piecewise constant, since every piecewise affine 
function can be approximated arbitrarily well in the L°° topology by piecewise constant func- 
tions and because the map associating to v(-) the curve t i— > M{T)dT is continuous with 
respect to the L°° topology (taken both in its domain and its codomain). 
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Assume that f (■) is constantly equal to w G R on [ti, ^2]- Since A2 — Ai, . . . , Aat — Aat-i are 
Q-linearly independent, then for every Sq G R the curve 

(so, 00) 3 ((Ai - A2)S, . . . , (Ai - AAr)s) 

projects onto a dense subset of the torus R^~-'^/27rZ^~-'^. Thus, there exist two sequences 

yjirn) y j^n^J ^(m) y g^^j^ ^^^^ 

(Ai - Xj)w^"'^ (mod 2%) — > (Ai - \j)w (mod 2%) for 2 < j < A^, 

(Ai - Xj)z^"^^ (mod 27r) — > (Ai - Aj)w (mod 27r) for 2 < j < n, 

(Ai - Aj>("*) (mod 27r) — ^ (Ai - Aj)m; + tt (mod 27r) forn + 1 < j < iV, 

as m tends to infinity. In particular the sequence of matrices e-«"^'"'^^^'i^(^)e'"^'"'^^^-' con- 
verges to e~"'"^''^'i?(^)e"''^''^' as m goes to infinity, while the sequence e~^*'"'^*'^'i?(^)e^''"''^^^'' 
converges, following the notations introduced in (4.7), to 

where we dropped the dependence on t of the different sub-matrices since v is constant on 

Fix 6 > 6. Consider a sequence {vk)keN in R+ (whose role will be clarified later) and define, 
for every A; G N, a finite increasing sequence {9f)i=o,...,k with 0^ = , Of^^ > 9f + 6(t2 — tij/k"^ 
for < Z < /c, and such that, for I > 0, 9i belongs to (ty^'^^)meN if I is odd and to {z^'^^)me]ss 
if / is even. Define, moreover, 

Tj =ti + ^-J^{t2 - h), 

for j = 1, . . . , /c + 1. 

Consider the continuous function Vk uniquely defined on [^1,^2] by the conditions 

Vk{Ti+'-^) = 9^ f0Tl = l,...,k, 

Vkit) = 6 if t G Uti(r, + ^2^,ri+i), 
iifc(t) = iiteuti{n,Ti + '-^). 

(See Fig. 1.) Define Vk as the piecewise constant function that coincides with 9'^ on [r^, Tj+i]. 
On each interval [Ti + {t2—ti)/k'^, Tj+i] the difference between Vk and Vk is bounded in absolute 
value by 6{t2 — ti)/k. Therefore, 



sup 



N) 



t G Ut, 



t2 — ti 



goes to zero as k goes to infinity. 

Since ||e~'^"^*^'-B*-^''e'^"^*^' || is uniformly bounded with respect to z/ G R and the measure 
of U^Tj, Tj + p-] goes to as /c goes to infinity, we have 

r (e-^(-)^'"^i?(^)e^^M^*"' - e-^^(-)^'"'i?(^)e*^(-)^'"')rfr '-^ uniformly on [t,, ^2] ■ 
Jti 
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Figure 1: The functions and 



Moreover, by definition of the sequences Q\ 



w 



and 



Qj^g Yias 



e -i).(r)A(^')^(7V)gS,(r)A(^-^^ ™ I M(r)rfr uniformly on [ti, tsj , 



and therefore e~'''='^*)'^^^''i?'^^)e"''''*^^^^'' converges in integral sense to M{t) on [ti,t2]- 

Finally, construct Uk as follows: for ti = define Uk on (^1,^2) as the derivative of Vk 
(defined almost everywhere), where the Wfc's correspond to the sequence of initial conditions 
Ufc = for every /c G N. Then, on the second interval on which v{-) is constant, use as a new 
set of initial conditions for the approximation procedure the values Vk = Vk{t2) and define 
again Uk as the derivative of Vk- Iterating the procedure on the finite set of intervals covering 
[0, T] on which v{-) is constant we obtain the required approximating sequence of piecewise 
constant control functions. □ 



4.4 Approximate controllability for the infinite-dimensional system 

Let Uk and Mk be defined as in Claim 4.3. The resolvent Rk(t, s) : of the linear 

time-varying equation 

Y = Mk{t)Y. 

converges, uniformly with respect to (t, s), to the resolvent s) : — > of 

Y = M{t)Y. 

(See, for instance, [4, Lemma 8.10].) Notice that R(t, s) preserves the norms of both the vector 
formed by the first n coordinates and the one formed by the last N — n. 

Let, for every k & 'N,ip^ be the solution of (4.2) corresponding to Uk- We have the following 
approximation property. 

Claim 4.4. For k large enough, 

||n„(e-'"=(^)^^^(T)) - n„(e-^(^)^^i)|| < 2e, (4.9) 
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where e is the positive constant which has been fixed at the beginning of Section 4-3. 

Proof. Define = e'^^^*)^ According to (2.3) (more precisely, its counterpart for 

equation (4.2)), the components qj{t) = e~^^^'"''^^^ (^l)^{t),(j)j'^ of q^{t) with respect to the basis 
of eigenvectors of A satisfy for almost every t E [0,T] 



(4.10) 



Therefore, the curves P''(t) = (g^(t), . . . , qn{t)Y and Q^{t) = (g^+i(t) 

P^(t) 



q%{t)Y satisfy 



Mk{t) 



P\t) 



H^{t) 
I\t) 



with ll-f^'^lloo < y/nfl (see (4.6)) and ||/*''||oo < C for C = C{N) large enough. 



Hence 



P\t) 



Pfc(t,o)n^(V^o) + / Rk{s,t) 



H\s) 

iHs) 



ds. 



Denote by 11^ the projection of C on its first n coordinates and let 



ds 



Since converges uniformly to R and the latter preserves the norm of the first n components, 
we know that, for k large, H-Z^'^Hoo < ^T.ynjl. Moreover, Rk{t, 0)Iliy{'ipo) converges uniformly to 

R{t,0)IlN{ipo)- In particular, according to the definition of R, {Rk{t,0)IlN{il>o)) converges 
uniformly to the solution of (On) corresponding to the control u and starting from Iln{ipo) = C,o- 
Since u drives system (S„) from ^o/||^o|| to Ci/||^i||, then it steers system (6„) from to 



-t)(T)A(") 



eidieoll/lleilD- Therefore, 



if k is large enough. Let us fix /j, small enough in order to have 

3Ty/nJi < e. 

Then, 

||n„(e-^'=(^)V(T))-n„(e-''(^)^V^i 



< 



n„(e-''^(^)^^^(T))-n„(e-''(^)^^i 



11^0 



ll^il 



116 

n„(e-''(^)^Vi) 



1161 



+ ||n„(e" 



^(m^.JII I - 1161 



116 



< 3Tv/7^+ I 11611 -1161 

< 25, 



provided that k is large enough. 



□ 
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As a consequence of Proposition 4.4, for k large enough the moduh of the first n components 
of ip^{T) are close to those of the first n components of ipi. The proposition below will be 
used to show that their phases can also be made as close as required by applying a suitable 
control on an arbitrarily small time interval. 

Proposition 4.5. For every £ > 0, every Ai,...,A„ G R, and every si G R there exist 
S2 > Si and w G R*^ with \\w\\ < e such that AjS2 = Wi mod 27r for every i = 1, . . . ,n. As 
a consequence, given a skew-adjoint discrete- spectrum control system {A,B, (0,5)), for every 
f 1 G R and every r > small enough there exists a constant control function u : [0, r] 
(5, +oo) such that every trajectory ip^-) of tp = uAip + BiIj corresponding to n(-) satisfies 
||n„('?/'(r)) — n„(e*'^'^'?/;(0))|| < e (where n„ denotes the orthogonal projection on the space 
spanned by the first n eigenvectors of A). 

Proof. The first part of the statement is a simple application of the Poincare recurrence 
theorem. Indeed, since the dynamics s i— > + s(Ai, . . . , A^) on the n-dimensional torus 
preserve volumes and distances, then the constant vector field (Ai,...,A„) is recurrent at 
every point of the torus, and in particular at the origin xq = 0. Therefore any neighborhood 
M of the origin is sent, after a suitably long time (which can be assumed to be larger than si), 
to another neighborhood of the origin isometric to M . Taking M equal to the ball of radius e 
centered at the origin, the first part of the claim is proven. 

In order to conclude the proof, fix a piecewise constant control function u : [0, r] (5, +oo) 
and a solution il) oi ip = uAip + Bip corresponding to u. Set 

q{t) = e- ^>^'^'^'^^{t) 

and notice that, according to (2.3), 15^(^)1 < for every j G N, where (pj denotes the j'-th 

eigenvector of A and qj = <^g, (pj^. Therefore, ||n„(g(r)) — n„(g(0))|| < Cr for some positive 
constant C independent of u and of tp. Then 

||n„(^(r)) - n„(e^^^^(0))|| = ||n„(e-''^^^(r)) - n„(^(o))i| 

= ||n„(e(^<'"W'^*^^O^g(r))-n„(g(0))|| 

< ||n„(e(^o^"W'^*-i)^g(^)) _ n„,(g(r))|| + Cr. 

Fix r < e/{2C) so that 

mn{i>{T)) - n.(e-^^(0))|| < ||n.(e(/o^'^W'^*-^O^g(r)) _ n.(g(r))|| + |. (4.11) 
Notice that 

n„(e(^o^"(*)'^*-"0^g(r)) = diag (^e'^i (/o^ "(*)'^*-''0 , . . . ^ e'^"(^o^"W'^*^"0 j n„(g(r)). 
The first part of the claim ensures the existence of V2 arbitrarily large such that if 

u{t)dt -vi = V2 (4.12) 
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then the norm of the matrix 



diag (^e*^i(/o""W'^*-''i), . . . , e*^"U'"W'^*^^'0 j - Id„ 

is smaller than e/2. Take V2 large enough to satisfy {vi + V2)/t > 6. Then u = (vi + V2)/t 
satisfies (4.12) and, because of (4.11), 

||n„(^(r))-n„(e^^^V;(0))||<e, 

independently of □ 
To conclude the proof of Theorem 2.4 we extend the interval of definition of the control 
function Uk introduced above by taking Uk{T + t)= u{t) for t G [0, r] where u is the control 
obtained by applying Proposition 4.5 with vi = v{T) — ffc(T), A = A, B = B, e = e. Then, 
for k large enough, the corresponding trajectory ip^ : [0,T + t] ^ S satisfies \\Iln{i''' {T + r)) — 
Uni'ipi)]] < Se and therefore, due to (4.5), 

||n„(V;'=(T + r))-^i|| <4£. 

It is now enough to show that ||'i^''"(r + r) — Iln{ip''(T + t))\\ can be made arbitrarily small by 
choosing a suitably small e. To this aim we notice that the inequality ||n„('?/;'^(T+r))|| > l—4e 
implies, for e < 1/4, that ||^/^*^(T + r) - n„(?A'=(T + r)) f < 1 - {1 - Aef = 8e - IGe^ and this 
concludes the proof of Theorem 2.4. • 



4.5 Lower bound on the steering time 

In this section we prove a lower bound on the steering time for a skew-adjoint discrete-spectrum 
control system without assuming that it satisfies the hypotheses of Theorem 2.4 nor any other 
controllability assumption. 

Proposition 4.6. Let {A,B,{0,6)) be a skew-adjoint discrete- spectrum control system. Fix 
ipo^ipi in S and e > 0. Then if a piecewise constant control u : [0,T„] — > (0,5) steers system 
(2.1) from ipQ to an e -neighborhood ofipi, then 

Tu > -^sup,,^ , (4.13) 

where (0fc)fcgN denotes the orthonormal basis of eigenvectors of A. 

Proof. Fix an initial condition ipQ in S, a piecewise constant control u : [0,T„] (0,5), 
and denote by : [0, T^j] — > Jif the corresponding solution of the system (2.1) satisfying 
7/;"(0) = ^Po. 

Write u as u{t) = I]J=o %Xfe,t,+i)(^) where = ti < ^2 < ■ ■ ■ < Wi = Tu and Mi, . . . ,m„ 
belong to (0,5). In the spirit of Section 4.1, associate to u the piecewise constant control 
u : [0, T^] R given by u{t) = YJj=Q '^jX[T„T,+i)it) with Uj = l/uj for all j = 1, . . . , n and re- 
defined by induction as ri = 0, Tj+i = Tj -\- {tj^i — tj)uj for j > 1. 

Define ip'^ : [0,Tjy] —>■ ^ as the solution of system (4.2) corresponding to v and satisfying 
?/;'^(0) = V'o- Define by = | , (pk) \ the modulus of the /c*^ coordinate of tp" . 
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By definition is absolutely continuous and equation (4.10) implies that 



oo / oo \ V2 

j=i \j=i J 

Applying the mean value theorem, one gets 

\\{r{0),4>k)\ - <TAB4>k\\. (4.14) 

Notice that T^, = X]j=i(^j+i ~ — (^n ~ ^i)^ — '^u^^ ^^^^ i^, 

Tu > \t,. (4.15) 



Since, by assumption, lip^iTi,) — ipi\ < e, then (4.14) implies 



Ty > SUPfcgN 



K0fc,^o)| - \ {^k,1pl)\\ 



Plugging this last inequality into (4.15), we obtain (4.13). □ 

We insist on the fact that this result is valid whenever system (2.1) is or is not approxi- 
mately controllable. 

Remark 4.7. When B is bounded, the same estimate as above is valid for other classes of 
controls (not only piecewise constant functions but also measurable bounded or locally inte- 
grable) as soon as we can define a unique solution of system (2.1) that satisfies (2.3). See 
Remark 2.2. 

Remark 4.8. It follows from (4.13) that, in general, approximate controllability does not imply 
finite-time approximate controllability. Indeed, if B(j)k tends to as /c goes to infinity, then 
for every T > the attainable set at time T from a given point i/jq is not dense in S since for 
every e G (0, 1), for k large enough, (pk is not e-approximately attainable from '0o in time T. 



5 Controllability for density matrices 
5.1 Physical motivations 

A density matrix (sometimes called density operator) is a non-negative, self-adjoint operator 
of trace class [32, Vol. I] on a Hilbert space. The trace of a density matrix is normalized to 
one. As a consequence of the definition a density matrix is a compact operator (hence with 
discrete spectrum) and can always be written as a weighted sum of projectors, 

oo 

P = ^Pjfjf*j^ (5.1) 
i=i 

where Pj G [0, 1], J2j Pj = ^j^* is the orthogonal projector on the space spanned by 

with </3j(-) = (v^j, ■)• Here {v^jjjeN is a set of normalized vectors not necessarily orthogonal. 
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The density matrix is used to describe the evolution of systems whose initial wave function 
is not known precisely, but only with a certain probability, or when one is dealing with an 
ensemble of identical systems that cannot be prepared precisely in the same state. More 
precisely (5.1) describes a system whose state is known to be ipj with probability Pj, j G N. 
Given an observable A (i.e. a self-adjoint operator, for instance the drift Hamiltonian) the 
mean value of A is Tr(py4) = Pj {Vj^ ^¥^j)^ where {ipj, Aipj) represents the mean value of 
the observable A in the state ipj. When for some /c G N we have = 1 and Pj = for every 
j 7^ k, one says that p describes a pure state, otherwise one says that p describes a mixed 
state. In the case of pure states, the physical description via the density matrix is equivalent 
to the one via the wave function. Notice that for a pure state Tr(p^) = 1 while for a mixed 
state one has Tr(p^) < 1. 

Without loss of generality it is possible to require that {y^jjjgN is an orthonormal basis 
(i.e. a basis of normalized eigenvectors of p). In this case {PjjjgN is the spectrum of p. 

The time evolution of the density matrix is determined by the evolutions of the states ipj, 
namely 

p(t) = U(t)p(0)U*(t) (5.2) 

where U(t) is the operator of temporal evolution (the resolvent) and U*(t) its adjoint. Notice 
that the spectrum of p{t) is constant along the motion. 

5.2 Statement of the result 

Fix 6 > and let {A, B, (0, 6)) be a skew-adjoint discrete-spectrum control system on a Hilbert 
space J^, ((^j)jgN an orthonormal basis of Jif (not necessarily of eigenvectors of A), {Pj}jQN 
a sequence of non- negative numbers such that Pj = 1, and denote by p the density 

matrix 

oo 

Definition 5.1. Two density matrices po and pi are said to be unitarily equivalent if there 
exists a unitary transformation U of such that pi = UpoU*. 

Obviously the controllability question for the evolution of the density matrix makes sense 
only for pairs (po, pi) of initial and final density matrices that are unitarily equivalent. Notice 
that this is a quite strong assumption, since it implies that the eigenvalues of po and pi 
are the same. Controllability results in the case of density matrices that are not unitarily 
equivalent have been obtained in the case of open systems (i.e. systems evolving under a 
suitable nonunitary evolution) in the finite-dimensional case. See for instance [9]. 

Next section is devoted to the proof of the following theorem. 

Theorem 5.2. Let po and pi be two unitarily equivalent density matrices. Then, under the 
hypotheses of Theorem 2.4, for every s > there exists a piecewise constant control steering 
the density matrix from po e- approximately to pi i.e. there exist /c G N, ti, . . . , > and 
ui, ...,UkE (0,5) such that setting V = e^kiA+UkB)^_ . _^^h{A+mB)^ ||p^_VpoV*|| < e, 

where || • || denotes the operator norm on M' . 
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Remark 5.3. As Theorem 3.4 is a particularization of Theorem 2.4 to the controlled Schrodinger 
equation, the hypotheses of Theorem 3.4 imply ^-approximate controllability of the corre- 
sponding density matrix. 

5.3 Proof of Theorem 5.2 

The proof uses the notations of Section 4. As noticed in Section 4.1, the theorem can be 
restated in terms of the evolution of the density matrix corresponding to the control system 
ip = {uA + B)^, u e (5, +oo). 

Fix po and pi unitarily equivalent and let U be such that pi = UpoU*. Write 

oo 

i=i 

with {Pj)j^N a sequence of non-negative numbers whose sum is one, and {(poj)j£N an or- 
thonormal basis of Jif. Then 

CXj 

i=i 

with (fij = U(^o,j for every j G N. 

Choose e > 0. Let m G N be such that 

The idea is to follow the strategy applied in the proof of Theorem 2.4 in order to simultaneously 
approximately steer m copies of system {A, B, (0, 6)) from (poj to ipij, j = 1, . . . , m. 

Let ?7 > be a small constant depending on m and e, to be fixed later. There exists 
n = n{7]) > m such that, for every j = 1, . . . ,m and for k = 0,1, 

\\<^k,j - n„v2fc,ill < v- 

By construction, when t] gets small, the two families {Iln^Pk,j)Y=i^ ^ ~ 0)-'-! tend to two 
orthonormal families. Hence, there exists a matrix M in SU (n) such that 

||M(n„(^o,i) -n„(^ijll (5.3) 

for J = 1, . . . , m provided that r] is small enough (and, consequently, n is large enough). 

Without loss of generality we may assume that B^^^ is connected. Claim 4.1 can be 
extended to the following result. 

Claim 5.4. The control system 

^ = (nAW + 5("))^, geU{n), (5.4) 

is controllable in the following sense: for any qq, qi in U{n), there exists a unitary complex 
number e*^ with < 6 < 271 /n, a time T > and a piecewise constant function u : [0,T] — 
(5, +oo) such that the solution : [0,T] U{n) of (5.4) with initial condition (7"(0) = go 
satisfies e^^g^{T) = gi. 
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Proof. Let us first assume tfiat at least one among and lias nonzero trace and hence 
does not belong to su{n). In this case the inclusion (4.3), with A^") = 21 and E^") = «B, 
implies that Lie(A("), E^'^)) = u(n). Classical controllability results for right invariant systems 
on compact Lie groups (see [24, 37]) ensure that the attainable set from go of (5.4) coincides 
with U (n) so that the claim holds with 6 = 0. 

It remains to consider the case in which the traces of A^"'"' and -B'-"^ are zero, i.e. A^""^ and 
S*^") belong to su{n). In this case (4.3) implies that Lie{A^^\ B^^^) = 5u(?t,), and therefore 
the attainable set from qq of (5.4) coincides with goSUln), the set of matrices of U{n) having 
the same determinant as go. Given a target gi there exists G [0, 27r] such that det{go) = 
e"""^ det(5'i) = det(e~*n (71). Hence the claim holds true with 6 = d/n. □ 

Let T > 0, n : [0, T] — (5, +00) and < 6' < 2TT/n be such that the control u steers system 
(5.4) from /„ to e*^M. Notice that, without loss of generality, 27r/?7, < e. 
Let /i > be a small constant to be fixed later. Fix G N such that 

\\{bjl)l>N\\l^ < 

for every j = 1, . . . , n. Let us apply Claim 4.3 to the control function u and denote by {uk)ke'N 
the sequence of piecewise constant control functions obtained in this way. Write, moreover, 
v{t) = j^u{T)dT and Vk{t) = j^Uk{T)dT. For every /c G N write Uk as 

Pk 

Uk{t) = J2 ^jX[t^,t^+,)ii)^ t e [0, T] 
i=i 

with = t\ < ■ ■ ■ < t^p^ = T and denote by the unitary transformation 
For every j = 1, . . . , m, 

l|n. (e(''(^)-'''=(^))^v,^o,.) - n„ {^,,) II < ||n„ (e(^(^)-^'=(^))^v,y,o.) - e^'M (n^v^o.) II 

+ ||e^^M (n„(^oj) - M (n„(^oj) II 

+||M (n„(^o,i) -n„((^i,j) II . 

The same computations as in Section 4.4 (cf. (4.9)) show that, for every j = 1, . . . ,m, 

||n„ (e(''(^)-''=(^))^V,^o.) - e^'M (n„^o.) || < 2e 
for /i small and k large enough. Since Q<9<2Tx/n<e, then, for every j = 1, . . . , m, 

Ije^^M {lin^o,j) - M {Tin^o,j) II < \e'' -l\<e. 
Hence, because of (5.3), for k large enough, for every j = 1, . . . , m, 

||n. (e(^(^)-^^(^))^V,v^o,.) - n„ (v^i,,) II < As. 
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Applying Proposition 4.5 we can, up to the extension of to a piecewise constant control 
defined on a larger interval, assume that 



for every j = 1, . . . , m. Therefore, 

oo 

llVfcpoVfc - Pill = W^Pj ((Vfe(^oj)(Vfc(^o,i)* - ^hMj) II 

m. 

^Pj- ((Vfcy?oj)(Vfc</?o,i)* - y'lj-V'tj) II + 2£ 



< 



< 



J2 (ll^t^ojll + ||^i,ill) ||Vfe(po,, - <^i,,|| + 26 
< 2(5£) + 2e = 125, 
provided that k is large enough. 

6 Examples 

6.1 Perturbation of the spectrum 

The scope of Section 6 is to show how the general controllability results obtained in the previ- 
ous sections can be applied in specific cases. In particular, we want to show how the conditions 
on the spectrum of the Schrodinger operator appearing in the hypotheses of Theorem 3.4 can 
be checked in practice. 

Let us adopt the notations of Section 3 for the domain fi, the wave function ip, and the 
uncontrolled and controlled potentials V and W . Throughout this section we assume that 
one of the hypotheses (i) or (ii) of Corollary 3.3 holds true. Thus, (A, S, U) is a well-defined 
controlled Schrodinger equation, where A = —i{—A + V) and B = —iW. 

The study of the examples below is based on the simple idea that, even if the hypotheses 
of Theorem 3.4 are not satisfied by the operators A and B, one can anyway ensure that they 
hold true for A^ = — i(— A + \^ + fiW) and -B^ = —iW for some /i in the interior of U. This 
is enough to conclude that the system ip = Aip + uBip^ u E U, is approximately controllable, 
since the replacement of {A, B) by (A^, S^) corresponds to a reparameterization of U that 
sends u into a new control u — ^ E U — fi and V into V + ^W. Although the spectrum of A^ is 
not in general exphcitly computable, we can nevertheless deduce some crucial properties about 
it by applying standard perturbation arguments. Theorem 6.1 recalls, in a simplified version 
suitable for our purposes, some classical perturbation results describing the dependence on /i 
of the spectrum of -A + V+fj.W. (See [25, Chapter VII, Remark 4.22], [33, §11.10, Theorem 1] 
and also [6].) 

Theorem 6.1. Let U he an open interval containing zero. Assume either that (i) Q is bounded, 
V,W belong to L°°(fi) or that (ii) fl = R'^, V belongs to Ll^i'R'^), W belongs to L°°(R^), 
lim||2.||^_i„oo ^(2;) = +00 and ini ^^-^d V (x) > —00. In both cases (i) and (ii) assume that each 
eigenvalue of the Schrodinger operator —A + V is simple. Denote by {\k)k&i the sequence of 
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eigenvalues of —A + V and by {4>k)k(^'N the corresponding eigenf unctions. Then, for any k in 
N, there exist two analytic curves : U C and $fc : f/ — > L^{Q) such that: 



• for any ix in U, (Afc(/i))fegN is the family of eigenvalues of A — V + jjW counted ac- 
cording to their multiplicities and ($fe(/^))feeN is an orthonormal basis of corresponding 
eigenf unctions; 



We check below that if the derivatives A'^(O) are Q-hnearly independent then for almost 
every fi & U the eigenvalues of —A + ^ + are Q-linearly independent. This fact is used 
in the following to apply Theorem 3.4 to situations in which the uncontrolled Schrodinger 
operator has a resonant spectrum. 

Recall that, in the notations of Section 4, for any pair of integers j, /c G N, 



In particular, A'^,(0) = W{x)\(j)k{x)\'^dx is equal to b^k- 

Proposition 6.2. Let U be an open interval containing zero and assume that VL, V and W 
satisfy one of the hypotheses (i) or (ii) of Theorem 6.1 and that the eigenvalues of —A + V are 
simple. If the elements of the sequence {bkk)keN (if^G Q-linearly independent, then for almost 
every fi in U the elements of {Ak{fi))ke'N ore Q-linearly independent. 

Proof. Let / G N and z = {zi, . . . ,zi) G Q'. Denote by the subset of elements fi in U 
such that Yl]=i ^j-^jif^) — 0- Since each /i i— Afc(/i) is an analytic function, then is either 
equal to U or to a countable subset of U. Since bu = A[{0), . . . ,bu = A;(0) are Q-hnearly 
independent, then = f/ if and only if 2; = 0. Hence, the union T = U;gN U^gQi ^^o ^2 has 
Lebesgue measure zero, since it is countable. By construction, if fi does not belong to T, the 
elements of (Afc(/i))fc6N are Q-linearly independent. □ 

The other crucial hypothesis of Theorem 3.4 is that fcjj+i 7^ for every j G N (or, 
more generally, that i?^"-* = {bjk)^j^^i is frequently connected, see Remark 4.2). By the same 
analyticity argument as above one checks that either such hypothesis is always false or it is 
true for almost every /i G f/. 

Corollary 6.3. Let U be an open interval containing zero and assume that Vt, V and W satisfy 
one of the hypotheses (i) or (ii) of Theorem 6.1 and that the eigenvalues of —A + V are simple. 
Assume moreover that the elements of the sequence (6fcfc)fceN are Q-linearly independent and 
that i?'-"'-* is frequently connected. Then the controlled Schrodinger equation associated with VL, 
V , W and U is approximately controllable for every U dU with nonempty interior. 



• Afc(O) = \k and <l>fc(0) = 0^; 




(6.1) 
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6.2 ID harmonic oscillator 



In this section we study the Schrodinger equation describing the evolution of the controlled 
one-dimensional harmonic oscillator, 

t^{t, x) = -^(t, x) + {x' - u{t)W{x)) iPit, x), (6.2) 

where ip is the wave function depending on the time t and on a space variable x G R = fi. 
Recall that u{-) is a piecewise-continuous function with values in a subset U of R. Notice 
that the potential corresponding to the uncontrolled Schrodinger operator is V{x) = x^. The 
control system (6.2) has been studied, among others, by Mirrahimi and Rouchon who proved 
its non-controllability in the case where W is the identity function (see [28]). 

As a consequence of Theorem 3.2, the spectrum of — A-l-^ is discrete. Its explicit expression 

is 

{Afc = 2A; + 1 I A; > 0} , 

and therefore A^+i — are Q-linearly dependent. Each A^ is a simple eigenvalue whose 
corresponding eigenfunction is 

where Hk{x) = (— l)'^e^^^e~^^ is the k^^ Hermite polynomial. 

In order to apply Corollary 6.3 we would like first of all to ensure that the elements 



Jkk 



f W{x)Hk{x)-^e-'''dx, k>0, (6.4) 



are Q-linearly independent. Notice that for W{x) = x (i.e., the non-controllable case pointed 
out by Mirrahimi and Rouchon), since each function 0^ is even, bkk = f Wcpl = 0. 

The existence of controlled potentials W for which the elements of (bkk)k£'N are Q-linearly 
independent can be easily inferred from the linear independence of the functions 0^. The 
proposition below provides some explicit W with such a property (and such that the corre- 
sponding Schrodinger equation is controllable). The potentials W will be chosen in L°°(R) 
and therefore, as already remarked in Section 3, the corresponding solutions in the sense (2.2) 
coincide with mild or strong solutions, depending on the regularity of the initial condition. 

Proposition 6.4. (1) IfW is even, then system (6.2) is not approximately controllable. (2) 
If W has the form W : x ^ e"^ +bx+c^ with a, 6, c G R such that a < and the two numbers 
a/1 — a and b are algebraically independent, then system (6.2) is approximately controllable, 
provided that U has nonempty interior. 

Proof. Since each function 0^ has the same parity as the integer /c, then 0^0^ has the same 
parity as the integer j + k. \{ W is even, then (6.1) shows that for every (j, /c) such that 
j + k is odd, bjk = 0. Applying Remark 4.2, one sees that the spaces spanned by the sets 
{0fc I k even} and {0^ | k odd} are invariant by the dynamics of system (6.2). In particular, 
there is no way to steer system (6.2) from 0i to a point 5-close to 02 if e is smaller than ^/2. 
This proves (1). 



21 



In order to prove (2) let us apply Corollary 6.3 (with U playing the role of U and R the 
role of U). Let W have the special form W : x ^ e""^ H-fox+c ^ multiplication of W by 



the strictly positive real number e^^"-^) we may assume without loss of generality that 
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4(a- 1) 

Using the specific expression (6.3) of (pk in the definition of hjk we can write 



(6.5) 



h,k = (-l)V,afc£e'^^'+^^+^i/fc(x)£-e-^'rfx, 



with CT; = 1/ y^/!2\7^, / = /c, j. Notice that Hk is a polynomial with rational coefficients and 
of degree /c, whose leading coefficient is equal to 2^. Integrating by parts j times, we get 

where Pj^k is a polynomial of degree j + k. Define {gi^Ym^Q through 

j+k 

m=0 

Each can be seen as the evaluation at 6 of a polynomial with coefficients in Q[a] 
whose degree is less than or equal to j. If m G {k, k + 1, . . . , k + j} then has exactly 
degree j + k — m and the coefficient corresponding to the monomial of order j + k — m is 

The renormalization of c performed above is such that {a—l)x'^+hx+c = (a— 1) (^x + 2^i„a j • 
Hence, the change of variables y = y/1 — a (^x + j yields 



■/T^Jr VvT^ 2(a-l) 

Due to the remarks made above on the coefficients of Pj^k, we have 

^ ' ^""^'^ -V'-^ + QUhy) (6.6) 



where Qj^k is a polynomial with coefficients in Q(a/1 — a) (d QH) and of degree smaller than 
j + k m its first variable. Notice that the coefficient multiplying in (6.6) is different from 
zero. 

For every m > the integral /j^e~^ y^dy is equal to zero if m is odd and to V (^^^) = 
2m'(m-)i a/tt if m is even, where F is the Euler gamma function. 
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Therefore, if j + k is even, 



bjk — ~^7=^=Sj,k{b) 
V 1 — « 



where Sj^k is a polynomial with coefficients in Q(a/1 — a) of degree exactly j + k. 

Since b is transcendental over Q(-\/l — a) then 7^ as soon as j and k have the 
same parity. Moreover, the elements of the sequence (A'^(0))fc>o = (bkk)k>o are Q- linearly 
independent. 

To conclude the proof let us check that each matrix {bjk)'^^^Q is connected. Fix j, k G 
{0, . . . , n}. We should prove the existence of a sequence ri, . . . , r; G {0, . . . , n} such that 
hjr^br^r2 ' ' ' b'n^^ribrik 7^ 0. If j and /c have the same parity then we are done since bj^ ^ 0. 
Otherwise, a simple computation and the normalization (6.5) show that 



7^0 



V2{1 - a)3/2 

and we can conclude by taking {ri,r2} = {0, 1}. □ 

6.3 3D potential well 

Consider the Schrodinger equation 

dijj 

i^U, x) = -Atp(t, x) + u(t)W(x)ip(t, x), (6.7) 
ot 

where the wave function depends on the time t and on three space variables xi,X2,X3 
with {xi,X2,X3) G (0,/i) X (0, /2) X (0, /a) = and satisfies the Dirichlet boundary condition 
'ipldfi = 0. Notice that the potential corresponding to the uncontrolled Schrodinger operator is 
V{x) = 0. For every W measurable bounded, solutions in the sense (2.2) coincide with mild 
or strong solutions, depending on the regularity of the initial condition. 
The spectrum of the Schrodinger operator is 



h.2 h.2 1.2 



Kmm = ^ it + it + tI \ki,k2,k3>l 



For the sake of simplicity, assume that {h^Y, (^1^3)^, and {hhY are Q-linearly independent, 
so that all the eigenvalues are simple and the perturbation result appearing in Theorem 6.1 
can be applied. (The case of multiple eigenvalues can be treated similarly, applying a refined 
perturbation argument as the one used in [6].) 

The normalized eigenfunction corresponding to \k^^k2M is given, up to sign, by 



^-3 

!>fci,fc2,fc3l3;i,a;2,X3' 



22 . fkiXin\ . {k2X2n\ . ( k^x^nX 
= sm — ; — sm — ; — sm — ; . 

^1/2/3 V ^1 / V ^2 ; V ^3 ; 



Vhkl 

Proposition 6.5. Let {hhY , {hhY , one? {hhY be Q-linearly independent and define W{xi, X2, x^) 
^axx^+a^xi+a^x-i yjj^ff^ Oi, 0^2, ^3 G R. Assume that ai, 02,0:3 are nonzero and that (vr/ai/i)^, 
(7r/a2^2)^; (7r/a3/3)^ are algebraically independent. Then the control system (6.7) is approxi- 
mately controllable. 
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Before starting the proof of Proposition 6.5 let us show the following technical result. 

Lemma 6.6. Let j3 be a real number transcendental over a field F with Q C F C R. Then 
the elements of the family ij^^j are F-linearly independent. 



Proof. Fix G N and distinct numbers qi,...,qN G Q \ {0}. Assume that for some 
fi,...fN in F 

^ 1 

We have to prove that /i = /2 = ■ ■ ■ = /at = 0. Multiplying (6.8) by n^j^(l + q^P) we get 

N /N-l \ 

=0 (6.9) 



fc=l \r=0 



where Skfi = 1 and, for r > 1, 



^ < jl < 32 < ■ ■ ■ < jr < N 
Jl , . . . , jr k 



By hypothesis, all coefficients of the left-hand side of (6.9), seen as a polynomial in /3, are 
equal to zero. Hence, (/i, . . . , fN)SN = (0, . . . , 0) where 

^ Sl,0 • ■ ■ Si,Ar_i ^ 

Sn = : : ■ 

\ SNfi ■ ■ ■ Sn,N-1 J 

A computation shows that det(S'Ar) = Ili<j<:k<N{(lk — Qj) ■ Hence, Sn is invertible and therefore 

/l = /2 = --- = /Ar = 0. " " □ 

Proof of Proposition 6.5. Theorem 6.1 and Fubini's theorem imply that the eigenvalues 
^kiMMif^) '^^ + '^^ ^ f*-*^ Dirichlet boundary value problem satisfy 



where 



axha2l2aM^^^kl + alll){An^kl + alll){An^kl + alll) 

^*'*^*'(lM(l*M(S*i+^ 

_ 64(e°i'i - l)(e"2'2 - l)(e":'':^ - 1)tt^ 
O — 



(ai/ia2/2a3^3)^ 

Let (3j = 47r^/(a|/|), j = 1, 2, 3. The Q-linear independence of the elements of (A'^^j..^^^^ (0)). 
is obtained from the expression above thanks to three nested applications of Lemma 6.6 with 
F = Q(/9i,/?2) and /? = /^g, F = Q(/5i) and /3 = /32, and F = Q and /3 = (3i. In order to 
complete the proof, let us check that every matrix S*^"^ is connected. (The conclusion then 
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follows from Corollary 6.3.) A straightforward computation shows that for every triples of 
positive integers (/ci, k^) and {hi, /12, h^) the integral 
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